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We measure the transmission matrix for microwave radiation propagating through ensembles of 
random waveguides of different lengths and scattering strength. We find that the probability dis- 
tribution of the transmittance T, the analog of the electronic conductance, evolves from Gaussian 
to log-normal and passes through a one-sided log-normal distribution in the crossover to Anderson 
localization. This asymmetrical distribution obtained for an ensemble with a value of the dimen- 
sionless conductance g of 0.37 is explained within the framework of an intuitive Coulomb gas model 
for the eigenvalues of the transmission matrix. 
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Mesoscopic fluctuations of conductance and transmis- 
sion limit the reproducibility of integrated circuits and 
the possibilities of imaging through opaque samples [HQ- 
At the same time, however, fluctuations of the field 
within disordered samples may also be exploited to lower 
the threshold for lasing [3| and sharpen focusing in ran- 
dom systems 0, H[. The importance of fluctuations in 
electronics was first recognized in conduction mediated 
by localized states [B-Q, but their impact was soon ob- 
served in universal fluctuations of conductance in samples 
in which electrons freely diffuse 13, II [ and for classical 
waves, such as light, microwave radiation and ultrasound 
[l2j |. Such fluctuations make it impossible to predict the 
characteristics of an individual disordered sample and, 
conversely, transport can not be characterized by the con- 
ductance or transmission in a single sample or even by 
averages over many samples. However, an account of 
fluctuations via probability distributions of conductance 
in ensembles of random samples with different physical 
dimensions could provide a full description of transport 
0-0]. The corresponding measurements have not been 
possible in electronic systems because of the difficulty 
of producing collections of statistically equivalent sam- 
ples. Measurements of fluctuations of conductance are 
typically made in single physical samples by varying the 
applied voltage or magnetic field (lp| . 

Though systematic measurements of the statistics of 
conductance over ensembles of random samples have not 
been made, extensive measurements of fluctuation of 
transmitted intensity and total transmission for a single 
incident channel 12[ and of correlation between channels 
[l3j have been carried out for classical waves. However, 
measurements of the distribution of the transmittance, T, 
which is the sum of total transmitted flux over all incident 
channels, have not yet been reported. T is equivalent to 
the conductance in units of the quantum of conductance 
and is often called the "optical" conductance since it is 
the sum of transmission coefficients that can in principle 
be measured optically [Hj|. The localization threshold 
lies at (T} = g = 1, where (...) indicates the ensemble 
average 



The fullest measurement of transmission is the field 
transmission matrix t whose elements tb a relate the field 
in outgoing channels b to the field in incident channels a, 
Eb = Yl n =i tb a E a lallll- Here Nis the number of prop- 
agating channels supported by the sample leads at a given 
energy or frequency. T may be expressed in terms of the 
elements of t and also in terms of the eigenvalues r„ of the 
matrix product ttf , T = J2a,b=i 1 2 = En=i T «- Opti- 
cal measurements of the transmission matrix have been 
exploited recently to focus light transmitted through 
strongly scattering media |5|. But the dimensions of the 
matrix were too small for the impact of mesoscopic cor- 
relation on transmission to be manifest. 

Simulations [l8| have found a highly asymmetrical dis- 
tribution of the logarithm of conductance just beyond the 
localization threshold i n ag reement with random matrix 
theory calculations lij 20] for quasi- ID samples. How- 
ever, because of the challenges of electronic and optical 
measurements, distributions of the electronic or optical 
conductance have not been made and the origin of the 
anomalous distribution which emerge in simulations and 
calculations have not been further clarified. 

In this Letter, we present measurements of the trans- 
mission matrix for microwave radiation propagation 
through ensembles of random waveguides with the val- 
ues of g ranging from 6.9 to 0.04. In the limits of diffu- 
sive and localized waves, the probability distribution of 
the optical conductance, P(T), is found to be a Gaus- 
sian and log-normal distribution, respectively. For the 
ensemble with value of g = 0.37 just beyond the localiza- 
tion threshold, a highly asymmetric distribution of In T 
is observed, which is log-normal for T below the peak 
of the distribution and exhibits an exponential tail for 
high values of T. This anomalous distribution for T can 
be explained by the repulsion between charges and their 
images associated with the transmission eigenvalues in a 
model developed in ref. [l7j|, based on Dyson's Coulomb 
model for eigenvalues of large random Hamiltonian ma- 
trix 2l|. For low values of g, distributions of T become 
log-normal and the variance of the InT approach -(InT) 
as predicted in the single parameter scaling theory (SPS) 
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of localization for one dimensional systems 0, 0| • 

The scattering media studied are collections of ran- 
domly positioned alumina spheres with a diameter of 
0.95 cm and refractive index 3.14 embedded in Styro- 
foam shells [22j]. Those samples are contained within a 
copper tube of diameter 7.3 cm giving an alumina sphere 
volume fraction of 0.068. The transmission matrix is ob- 
tained by measuring tb a between arrays of points on the 
incident and output surface for microwave radiation po- 
larized along the orientation of the wire antenna with 
use of a vector network analyzer. The dimensionality of 
the measured transmission matrix matches the number 
of propagating waveguide modes supported by the cop- 
per tube, with N ~ 66 in the frequency range 14.7 to 
14.94 GHz and N ~ 30 in the frequency range 10-10.24 
GHz. The wave is localized in the lower and diffusive in 
the upper frequency range. The sample tube is rotated 
and vibrated after each measurement of the full trans- 
mission matrix to produce a new realization of the ran- 
dom sample. Measurements of spectra of the transmis- 
sion matrix are made for sample lengths L — 23,40,61 
and 102 cm. The impact of absorption on the statistics 
of transmission is removed by Fourier transforming the 
field spectrum into the time domain and multiplying the 
time signal by exp*/ 2Ta , where t is the time delay and r a 
is the absorption rate [221 ] . 
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FIG. 1. (Color online) Probability distribution of optical 
conductance for two random ensembles with values of g — 0.37 
and 0.04, respectively, (a) P(lnT) for g = 0.37 (red dots) and 
g = 0.04 (green asterisk). The solid black line is a Gaussian 
fit to the data, where for g = 0.04, all the data points are 
included and for g = 0.37, only the data to the left of the 
peak are used in the fit. (b) P(T) for g=0.37 in a semi-log 
plot with an exponentially decaying tail. 



Measurements of the probability distribution of InT, 
P(lnT), for random ensembles with values of the dimen- 
sionless conductance of g — 0.37 and 0.04 are presented 
in Fig. la. P(lnT) is seen to be Gaussian for g = 0.04 
and highly asymmetrical for g = 0.37. The low transmis- 
sion side of P(lnT) is well fit by a Gaussian distribution, 
and to fall sharply above the peak. For T > 1.1, P{T) is 
seen in Fig. lb to fall exponentially in accord with ref. 
201. 




FIG. 2. (Color online) Charge model of transmission eigen- 
values and conductance, (a) Average positions of charges and 
their images with respect to different positions of the first 
charge xi in the random ensemble with g — 0.37. The dashed 
lines show the average positions of the charges for this en- 
semble, (b) Average positions of charges vs. InT in the same 
ensemble. The curly brackets indicates the averaging is over 
subset of transmission matrices with the specified value of 
InT. 

These results can be understood with the aid of the 
charge model in which transmission eigenvalues r„ are as- 
sociated with positions of parallel line charges at x n and 
their images at — x n embedded in a compensating con- 
tinuous charge distribution [ItJ • The transmission eigen- 
values are expressed as, r„ = 1 / cosh 2 (x„ ) . The repul- 
sion between two parallel lines of charges with same sign 
with potential In \xi —Xj\ mimics the interaction between 
eigenvalues of the random matrix, while the oppositely 
charged jellium background provides an overall attrac- 
tive potential that holds the entire structure together. 
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FIG. 3. (Color online) Changing of distribution of optical conductance with respect to different values of g and the underlying 
transmission eigenvalues, (a) P(T) for three values of g = 6.9, 0.61 and 0.37. A Gaussian fit to P(T) for diffusive waves with 
g = 6.9 is shown, (b-d) A typical spectra of transmission eigenvalues r n and T for (b) g = 0.37, (c) g = 0.61 and (d) g = 6.9. 



The average spacing between neighboring lines of charge 
is Z//£ while the background charge provides a screen- 
ing length of \JL/^, where L and £ are the sample and 
localization lengths, respectively. The spacing between 
charges is smaller than the screening length for diffusive 
waves, for which L < £, so that the Coulomb interaction 
is not effectively screened and the charges experience a 
logarithmic repulsion. This is the origin of universal con- 
ductance fluctuations for diffusive waves 111. The re- 



pulsion between the first charge x\ associated with the 
highest transmission eigenvalues t\ and its image placed 
at -x\ naturally provides a ceiling for t\ of unity. 

In Fig. 2a, we plot the variation of the average posi- 
tions of the charges for different values of the position of 
the first charge x\ in the random ensemble with g = 0.37. 
The average spacing between x\ and x-i increases as the 
value of X\ decreases. At the same time, the spacing be- 
tween X2 and X3 hardly changes as x\ moves towards the 
origin and their average position is little changed. This 
reflects the tendency to heal large fluctuations in charge 
positions for more remote charges and is a manifestation 
of electrical screening in this model. 

The source of the sharp cutoff in P(lnT) can be seen in 
the context of the charge model by examining the spacing 
of charges for different values of In T shown in Fig. 2b to- 
gether with the average positions of charges for different 
values of x±, shown in Fig. 2a. A relatively high value 



of T can only be obtained when the first charge is near 
the origin. This is an unlikely circumstance because this 
charge is strongly repelled by its image. P(T) would be 
expected to fall especially rapidly for values of T above 
unity since this would requires two charges along with 
their images to be close to the origin. 

In Fig. 3, we present P(T) for g = 0.37,0.61,6.9 as 
well as a spectrum of transmission eigenvalues for a single 
realization for each ensemble. For diffusive waves with 
g = 6.9, many transmission eigenvalues contribute to T 
leading to a Gaussian distribution for P(T), as seen in 
Fig. 3a. For g ~ 1, we have seen that transmission 
is dominated by the highest transmission eigenvalue, t\ 
except when T is comparable or greater than 1. 

Deep in the localized regime, g <C 1, T is essentially 
given by T\. In this limit, T ~ t\ ~ 4cxp(— 2x\) and x\ 
is hypothesized to follow a Gaussian distribution. This 
leads to the log-normal distribution for T seen in Fig. la 
for g = 0.04, as predicted by SPS theory. SPS further 
predicts that in the limit of small g, var(lnT) = a 1 = 
— (lnT), so that P(lnT) depends only upon the single 
parameter -{\nT)/L = l/£. The ratio of a 2 and -(lnT) 
vs £./£ is plotted in Fig. 4 and seen to approach unity 
for L » ( ~ 24cm in accord with SPS. 

In conclusion, we have observed the joint distribution 
of T and the transmission eigenvalues r n in the Anderson 
localization transition. We have demonstrated that the 
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FIG. 4. (Color online) Hypothesis of single parameter scaling 
theory of localizaion. The ratio of variance of InT and -(InT) 
vs. the ratio of sample length L to the localization length £. 
The dashed green line is the prediction from SPS for L S> £. 



Coulomb charge model including image charges provides 
an intuitive basis for understanding the statistics of con- 
ductance. Since the statistics of transmission for a single 
incident channel, Tb a and T a can be measured for those 
samples, the full range of statistic of waves of nonintcr- 
acting particles can be studied in a single system. This 
demonstrates the power of a unified approach to meso- 
scopic physics encompassing electronic conductance and 
statistical optics. 
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